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In the real-time thermal field theory, the standard expression of shear viscosity for the nucleonic 
constituents is derived from the two point function of nucleonic viscous stress tensors at finite 
temperature and density. The finite thermal width or Landau damping is traditionally included in 
the nucleon propagators. This thermal width is calculated from the in-medium self-energy of nucleon 
for different possible pion-baryon loops. The dynamical part of nucleon-pion-baryon interactions are 
taken care by the effective Lagrangian densities of standard hadronic model. The shear viscosity to 
entropy density ratio of nucleonic component decreases with the temperature and increases with the 
nucleon chemical potential. However, adding the contribution of pionic component, total viscosity 
to entropy density ratio also reduces with the nucleon chemical potential when the mixing effect 
between pion and nucleon components in the mixed gas is considered. Within the hadronic domain, 
viscosity to entropy density ratio of the nuclear matter is gradually reducing as temperature and 
nucleon chemical potential are growing up and therefore the nuclear matter is approaching toward 
the (nearly) perfect fluid nature. 


I. INTRODUCTION 


The recent hydrodynamical [l|,[3 as well as some trans¬ 
port studies ii have indicated about an (nearly) ideal 
fluid nature of nuclear matter, which may be produced 
in the experiments of heavy ion collisions (HIC) like Rel¬ 
ativistic Heavy Ion Collider (RHIC) at BNL. The hydro- 
dynamical calculations became very successful in explain¬ 
ing the elliptical flow parameter, V2 from RHIC data P- 
Q only when they assumed a very small ratio of shear 
viscosity to entropy density (77/s) for the expanding nu¬ 
clear matter. When some recent studies @t 11| (see also 
Ref. 0) show that rj/s may reach a minimum in the 
vicinity of a phase transition, then some special atten¬ 
tions are drawn to the smallness of this minimum value 
with respect to its lower bound {rj/s = ^), commonly 
known as the KSS bound [l^. In this context, the tem¬ 
perature (T) dependence of 77/5 is taken into account in 
some recent hydrodynamical calculations [IM3 instead 
of its constant value during the entire evolution. Niemi 
et al. [I3 have interestingly observed that the V2{pt) of 
RHIC data is highly sensitive to the temperature depen¬ 
dent 77/s in hadronic matter and almost independent of 
the viscosity in QGP phase. This work gives an addi¬ 
tional boost to the microscopic calculations of 77/s of the 
hadronic matter in the recent years (Tsl - I^ . though his¬ 
torically these investigations are slightly old (35l - l40l| . 


Except a few [2J, , most of the microscopic cal¬ 

culations are done in zero baryon or nucleon chemical po¬ 
tential {fiN = 0). Along with the T dependence of 77 or 
77/s, their dependence on the baryon chemical potential 
should also be understood in view of the future experi¬ 
ments such as FAIR. In the work of Itakura et al. 
and Denicol et al. [3l|, we notice that the 77/s is reduced 
at finite baryon chemical potential, whereas Gorenstein 
et al. [3^ observed an increasing nature of 77/s with 
Itakura et al. have obtained 77 by solving the relativistic 


quantum Boltzmann equation, where phenomenological 
amplitudes of hadrons are used in the collision terms. 
Denicol et al. have calculated the 77 at finite T and pN 
by applying Chapman-Enskog theory in Hadron Reso¬ 
nance Gas (HRG) model, whereas Gorenstein et al. have 
taken a simplified ansatz of 77 to estimate 77/s in the van 
der Waals excluded volume HRG model. Similar to the 
ansatz of r]{T) taken by Gorenstein et al, the 77 itself 
increases with increasing temperature in Ref. [24| . but 
their 77/s are exhibiting completely opposite nature of 
T dependence. Therefore, the behavior of the 77/s may 
largely be influenced by the T dependence of entropy 
density s. 


Motivating by these delicate issues of shear viscosity at 
finite /ijv, the present manuscript is concentrated on the 
matter with nucleon degrees of freedom at finite T and 
fj-N- The nucleons in the medium can slightly become 
off-equilibrium because of their thermal width or Lan¬ 
dau damping, which can be originated from the nucleon 
thermal fluctuations into different baryons and pion. The 
inverse of nucleon thermal width measures the relaxation 
time of nucleon in the matter from which one can esti¬ 
mate its corresponding shear viscosity contribution. 


In the next section, the one-loop expression of 77 for 
nucleon degrees of freedom is derived from the Kubo re¬ 
lation, where a finite thermal width is traditionally in¬ 
cluded in the nucleon propagators. This standard ex¬ 
pression of 77 can also be deduced from relaxation time 
approximation of kinematic theory approach. In the real¬ 
time thermal field theory, the nucleon thermal width from 
the different pion-baryon loops are calculated in Sec. 3, 
where their interactions are determined from the effec¬ 
tive hadronic model. In Sec. 4, the numerical results are 
discussed followed by summary and conclusions in Sec. 5. 
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KUBO RELATION FOR SHEAR VISCOSITY 
OF NUCLEAR MATTER 


II. 


From the simple derivation of Kubo formula [42 
let us start with the expression of shear visc osity for nu¬ 
cleonic constituents in momentum space (^ . l4f 


1 ^r,(90,^) 

riM = — lim — - - 

20 qQ,q-^0 Qq 


( 1 ) 


Hence, Eq. o can broadly be redefined as 


1 Imn«(go,<7) 

= Td - 

10 qo,q^o qo 


J_ e(go)Imn(go,^ 

10 go ,?—>0 


90 


1 tanh(,5go/2)Imn“(go,^) 

— lim - 

10 qo,q->-0 go 


(9) 


Using ([5]) in the 11 component of (O and then applying 
the Wick’s contraction technique, we have 


where 

Ariiqo,^ = J ( 2 ) 

is the spectral representation of two point function for 
nucleonic viscous-stress tensor, and 


nii(9) 


where 


I d*xe^‘^^{Ti; {xMi;{x)mrd'^m )p 

i J -^^N{q,k)Dii{k)Dii{p = q + k) , ( 10 ) 


p-PHA 

{0)p = Tr- with Z = Tre"^^ (3) 

Z 

is denoting the thermodynamical ensemble average. The 
energy momentum tensor of free nucleon is 

Tpa = + 

= -gpaC -I- ii^jpdai’ , (4) 

and hence the viscous stress tensor will be 


N{q,k) = -lNdZTr['y^{q+k)''{^+^+mN)^pka{^+mN)] ■ 

( 11 ) 

This self-energy function, nii(g) for NN loop can dia- 
grammatically be represented by Fig. [TJA). In the co¬ 
moving frame, i.e., for u = (1, 0), the N{q, k) becomes 


N{q,k) = -Im 


32 ^ 

y{fco( 9 o + ko)}{k ■ (g -h k)} 


-4 


{k- (q + k)} 


2 , P(g-yfc )2 


.( 12 ) 


— tpijTpa 

= tpZ'i'fPlpdcri’ ( since tfZgpa = 0 ) , (5) 

where 

= gf^''- uV . ( 6 ) 

In real-time formalism of thermal field theory, the ensem¬ 
ble average of any two point function always becomes a 
2x2 matrix structure. Hence, for viscous-stress tensor, 
the matrix structure of two point function becomes 

Uabiq) = i J d^xe^‘>^T,7rp,{x)7r^^’^{0))f , (7) 

where the superscripts a, b{a, b = 1,2) denote the thermal 
indices of the matrix and Tc denotes time ordering with 
respect to a symmetrical contour 1 ^ in the complex 
time plane. 

The matrix can be diagonalized in terms of a single 
analytic function, which can also be related with the re¬ 
tarded two point function of viscous-stress tensor. The 
retarded function n^(g), diagonal element n(g) and the 
spectral function H,j(g) are simply related to any one of 
the components of nab(g). Their relations with II com¬ 
ponent is given below 

Ar,{q) = 2Imn^(g) = 2e(go)Imn(g) 

= 2tanh(^)Imnii(g) . (8) 


In the above equations, Jjv = 2 is the isospin degeneracy 
of nucleon. 

In Eq. pOL is scalar part of 11 component of the 
nucleon propagator at finite temperature and density. Its 
form is 

D“(fc) = —- 1——- - 27:iFkiko)S{P - m%) 

k^ — + ig 

with Fk{ko) = n^0{ko) +n'j:9{-ko) 

1 . 1 - , < 

2io^ ^ ko — + ig kg — uj^ — ig 

_ . ^3^ 

ko+uj^ — ig ko+uj^ + ig’ 

where n^(w^) = _|_ ]^| jg Permi-Dirac dis¬ 

tribution function for energy uj^ = ^JkP-+ . Here 
the ± signs in the superscript of Uk stand for nucleon 
and anti-nucleon respectively. Among the four terms 
in Eq. m, the first and the second terms are associ¬ 
ated with the nucleon propagation above the Fermi sea 
and the propagation of its hole in the Fermi sea respec¬ 
tively, while the third and fourth terms represent the 
corresponding situations for anti-nucleon. The full rela¬ 
tivistic nucleon propagator, thus, treats the particle and 
anti-particle on an equal footing and all possible singu¬ 
larities (nucleon, hole of the nucleon, anti-nucleon and 
hole of the anti-nucleon) are automatically included. 
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N (k) n (I) 



(A) (B) 

FIG. 1: Diagrammatic representation of NN loop is shown 
in (A), where double lines stand for effective N propagators, 
which contain their thermal widths F. The diagrammatic 
representation of nucleon self-energy for tvB loop is shown in 
(B) from where F can be determined. 


After doing the ko integration of Eq. dTU)) and then 
using it in Eq. (ED, we have 


Vn 


1 f dPk (—ttN) 

(2^)3 


- - - S{qo-UJk 

qo 


,N\ 


- - - Hqo + Wfc - ) -b .. 


qo 


(14) 


where N = N{ko = ±uj^,k,q) and ujp = 

\/iq + kY + m%. 

The two ^-functions will be responsible for generating 
the Landau cuts {—q < qo < q), where the Imn'^(g) will 
be non-zero. However, there will be two more ^-functions 
(not written explicitly), which are not important for the 
limiting point qo,q —» 0 since th ey wil l generate unitary 
cuts (—oo < qo < —y/(p + 4m^ and \/(p + 4to^ < qo < 
oo). 

Using the identity 


— TrS{x) = Im 
in Eq. da, we have 


r]N = — lim Im 

10 (jo.g —>0 


1 


lim 
rjv-s-o X + iT 


N 


(Pk N 
( 2^)3 4 wf r ™o 


l-^fc (^fc ) + ^p (-go +^fc )}/go 
{qo - wf + uj^) + iT]s, 

{^t {^k ) - (go +^k)}/qo \ 

{qo + — ojp) + i^N J 


(15) 


(16) 


We will continue our further calculation for finite value 
of Ejv to get a non-divergent contribution of 77 at. In¬ 
cluding thermal width E n for constituent particles (here 
nucleons) of the medium is a very well established tech¬ 
nique [Ifl [ 23,113 in Kubo approach to remove the di¬ 
vergence of ryAT as well as to incorporate the interaction 


scenario, which is very essential for a dissipative system. 
The interaction scenario is coming into the picture by 
transforming the delta functions to the spectral functions 
with finite thermal width. The thermal width (or col¬ 
lision rate) Eat of the constituent particles reciprocally 
measures the shear viscosity coefficient, which is approx¬ 
imately equivalent to the quasi particle description. 

In the limiting case of go, 9 0, we get —>■ uj^ and 

therefore Eq. dMD is transformed to 



d^k {-No) 
(27r)3 4wf ^Eat 


[h + h] 


where 


No = lim N{ko = ±w^, k, q) 

qo,q^0 


and 


/ 2,3 = lim 

90->0 


{uJk ) ± nj{Tqo +UJ^)} 
qo 


(17) 


(18) 


(19) 


In the above Eq. dUD, one can notice that the limiting 
value of / 2,3 is of the 0/0 form. Therefore, we can apply 
the L’Hospital’s rule, i.e., 




^{TnJ(a;f)±nJ(Tgo+a.f)} 

lim ^- 

I3[nl{l-r4)\ , ( 20 ) 


since 


— {±nj(wg = =Fgo + )} 


Yivcl —{±n^{ujq = ipgo +wf)} 

90—>0 dqo ^ 


P dqo ^ _ 

|g/3(iv,±/Jiv) _|_ ]^j.2 

— (ip) ) 

{g/5Ur±/^v) _|_ ip 
^[<(l-^fc)] -( 21 ) 


Again, in the limiting value of go, t 0, Eq. (fT^ can be 
simplified to 


= -In 


16k^ 

3 


( 22 ) 


Hence, using the above results, the Eq. dnD becomes 


Vn = 


SpiN f d?k k'^ 
15 

PIn 


(27r)3 4u;f Ejv 
k^dk 


K(i-^fe) + <(i- 0 ] 


157r2 




[%(l-’^fc)+’^fc(l-’^fc)] ■ (23) 


N 


This is the one-loop expression of shear viscosity for the 
matter with nucleon degrees of freedom in the Kubo ap¬ 
proach. Though there are possibility of infinite number of 
ladder-type diagrams, which are supposed to be of same 
order of magnitude (©(l/E^v)) like the one-loop, they 
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Baryons 

Jb 

Ib 

Etot 

Vg^NiT (B.R.) 

f/m^ 

A(l232) 

3 + 
2 

3/2 

0.117 

0.117 (100%) 

15.7 

A* (1440) 

1 + 

2 

1/2 

0.300 

0.195 (65%) 

2.5 

A* (1520) 

3 - 

2 

1/2 

0.115 

0.069 (60%) 

11.6 

A* (1535) 

1 - 
2 

1/2 

0.150 

0.068 (45%) 

1.14 

A* (1600) 

3 + 

2 

3/2 

0.320 

0.054 (17%) 

3.4 

A*(1620) 

1 - 
2 

3/2 

0.140 

0.035 (25%) 

1.22 

A* (1650) 

1 - 
2 

1/2 

0.150 

0.105 (70%) 

1.14 

A* (1700) 

3 - 

2 

3/2 

0.300 

0.045 (15%) 

9.5 

A* (1700) 

3 - 

2 

1/2 

0.100 

0.012 (12%) 

2.8 

A*(1710) 

1 + 
2 

1/2 

0.100 

0.012 (12%) 

0.35 

A* (1720) 

3 + 
2 

1/2 

0.250 

0.028 (11%) 

1.18 


TABLE I: From the left to right columns, the table contain the 
baryons, their spin-parity quantum numbers Jg, isospin Ig, 
total decay width Ttot, decay width in Ntt channels T g^N-w or 
r g{mg) in Eq. (1321) (brackets displaying its Branching Ratio) 
and at the last coupling constants jIm-^. 


will be highly suppressed . As we increase the number 
of loops, the number of extra thermal distribution func¬ 
tions will also appear in the shear viscosity expression 
and hence their numerical suppression will successively 
grow. On this basis, the one-loop results may be consid¬ 
ered as a leading order results. One can derive exactly 
same expression from relaxation time approximation in 
kinetic theory approach. 


III. CALCULATION OF NUCLEON THERMAL 
WIDTH 

Now, our next aim is to calculate the thermal width of 
nucleon r^r, which can be estimated from the retarded 
component of nucleon self-energy (E^) at finite temper¬ 
ature and density. Their relation is given by 


B stand for different higher mass baryons including nu¬ 
cleon itself. In this work, all possible 4-star baryon res¬ 
onances with spin 1/2 and 3/2 are considered. These 
are A(980), A(1232), Ai*(1440), A*(1520), A*(1535), 
A*(1600), A*(1620), fV*(1650), A*(1700), A*(1700), 
A*(1710) and A*(1720), where masses (in MeV) of the 
baryons are given inside the brackets. The nucleon self¬ 
energy for 7ri? loop is shown in diagram [IJB) and its 11 
component can be expressed as 

E^^(fc,r,/XAf) = -*y 

Dii{u = k - I, niB, T, ^xn) , (25) 

where Dii{l, m,r, T), Dii{u = k—l,mB,T,^N) are scalar 
part of pion and baryon propagators at finite temperature 
and density. The L(k, 1) contains vertices and numerator 
parts of the propagators. The chemical potential of all 
baryons are assumed to be the same as nucleon chemical 
potential Similar to Eq. ([81), this 11 component is 
also related with retarded component as 

ImE^(/c) = coth ImEii(fc) . (26) 


Performing the Iq integration in (12511 and then using the 
relation (051) , we get the imaginary part of retarded self¬ 
energy. 


ImE^(fc) = TT 


d^l 


:[L{lo = u;iJ,k) 


(27r)3 4wfw^ 

[{1 -b nz(wn - n+{ko - uif)}S{kQ - w/ - wf) 
+{-ni{ujf) - n-{-ko + ujf)}6{ko - uf + wf)] 
+L{lo = -w/,/,fe)[{ni(u;n 
+n^ (fco + uji)}S{ko + iOi — 10^) + {—1 
-ni{uf) + n-{-ko - ujf)}S{ko + ujf + wf)]] , 

(27) 


where u!^ = y {k — l)'^ + m^, and n; are respectively 
Fermi-Dirac and Bose-Einstein distribution functions. 
The regions of different branch cuts in kQ-axis are (—oo 

to — \Jk?' + {rriT^ + ttib)^ ) for unitary cut in negative kg- 

axis, (—+ (ruB — m-nY to ^Jk‘^ + (tob — ) for 

Landau cut and (-b (m^ -I- tob)^ to oo ) for unitary 
cut in positive feg-axis. These are representing the dif¬ 
ferent kinematic regions where the imaginary part of the 
nucleon self-energy becomes non-zero because of the dif¬ 
ferent 5 functions in Eq. (1271) . The E m for all ttB loops 
(except the 'kN) are coming from the Landau cut contri¬ 
bution associated with the third term of Eq. (EH), which 
can be simplified as 


r^(fc, T, ^XM) = -ImE«(fco = , fc, T, ^MN) ■ (24) 

Eat = 

During the propagation in the hot and dense nuclear mat¬ 
ter, nucleon may pass through different ttB loops, where L(Zo = —w, I = \/lo^ — m^, ko = , k) ,(28) 
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where ni{w) = — 1 }, n+(w^ + w) = 


l/{g/3(c^+a;j, -Miv) 1}^ ^ ^ 

4m‘im'i, 1 7^)2 _ ^2 


W = \ \ - 


4 

o 2 ^2 


2mj, ± with 


j ^4 ^ and = mj^ + — m'^. 

The effective Lagrangian densities for BNtt interac¬ 


tions are given below 
f _ f 4^,5 

C- = —i’Br 


r f -fk' 
C = —Ipg 



interaction Lagrangian densities and second is to include 
the effect of these baryon widths (Ts) on the nucleon 
thermal width Tn- Using the Lagrangian densities, the 
vacuum decay width of baryons B for Ntt channel can 
be obtained as 


I - 

ipNdi^TT + h.c. for Jg = - 


ipNdf^n -I- h.c. for Jg = - , (29) 


T^BiruB) = 


T^BiruB) = 


2Jb + 1 


f 


\P' 


— ^^[2mB\Pc 


2TTmB 


- PrriN)] for Jb = 


V 


Ia 


,A*—^7tN 


TO, 


where coupling constants / /TOtt for different baryons have 
been fixed from their experimental vacuum widths in Ntt 
channel. With the help of the above Lagrangian densi¬ 
ties, one can easily find 




2JB-bl 
^ Ptoat] for 


^ IPemP 

SirruB 


Ju = 


(32) 


/ 


L{k,l) = - 



Vto^ 




Vto^ 

l 

2 

1 


3m?g 

-b 

1 

-( 


for Jg = - 


1 


3m B 


[k-iY[k-iY 


where \pcm\ = - - - - 2 mB - 

wYm = \/|PcmP + rnff. The isospin factors are 
= 3 and Ia,a*^-wN = 1 for the Ntt decay chan¬ 
nels of N* and A* (or A) respectively. 

Now, the rAT in Eq. (l28l) can be convoluted (see e.g. 
Refs, [li,!!^) as 


YYk-iY-{k-iYin> for j^ = 


r,{,ns) = YL 

Nb = 


(30) 


dMBAB{MB)T n{Mb) , 

pmB+2T B{mB) 

I Ab{Mb) , (33) 

tub— 2 r B^rnB) 


For simplification the coefficients of 7 ° and 1 are taken 
as in Ref. [i^ and their addition gives 


where 


Ab{Mb) = -Im 

TT 


Mb — mB + iTB{MB)l2 


(34) 


Hk,l) = - 


f 




— Pm^mB} for 


Jb = ^ 


L{k,l) = - 


TO, 




3to^ 


-to^to|} (fco - ^0 + PrriB) for Jb = ^ 


(31) 


is vacuum spectral function of baryons for their vacuum 
decay width in Ntt channel. Replacing baryon mass mB 
by its invariant mass Mb in Eq. (1321) , one can get the off- 
mass shell expression of T b{M). The values of coupling 
constants //to^, which are fixed from the experimental 
values of baryon decay width in Ntt channels [5l|, are 
shown in a Table (I). 


IV. RESULTS AND DISCUSSION 


The isospin part of the Lagrangian densities are not writ¬ 
ten in the Eq. (EH). The isospin structure for \ 

and JY = should be ijjf ■ ^ijj and ipf ■ Trip respec¬ 
tively, where T is the spin 3/2 transition operator and 
f is the Pauli operator. This issue is managed by mul¬ 
tiplying appropriate isospin factors with the expressions 
of corresponding loop diagrams. The isospin factor for 
ttN or ttN* loop is /Ar->TrAf,Af* = 3, whereas for the ttA 
or ttA* loop, Iat^-ti-a.a* = 2. 

All baryon resonances have finite vacuum width in Vtt 
decay channel. The calculations of these decay widths are 
very essential in the present work for two reasons. First 
is to fix the coupling constants f/m^r for different BNtt 


Let us first take a glance at the invariant mass dis¬ 
tribution of imaginary part of nucleon self-energy for 
different ttB loops. Fig. ED shows the results for 
baryons B = 7V(940), A(1232) (upper panel) and B = 
A*(1620), V*(1650), V*(1720) (lower panel), whereas 
Fig. © displays the results for baryons B = V*(1440), 
V*(1520), A*(1600) (upper panel) and B = V*(I535), 
A*(1700) (lower panel). The numerical strengths for 
B = V*(I700) and V*(I7I0) are too low to display with 
the other baryons. These results are obtained by replac¬ 
ing ujY = -I- m% by LUk = Vk'^ + M^ in Eq. (l28l) 

(dashed line) and (155)) (solid line) for the fixed values of 
k = 0, pn = 0 and T = 0.150 GeV. From the sharp 
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FIG. 2: Imaginary part of nucleon self-energy for different 
■kB loops are individually shown before (dashed line) and 
after (solid line) foldiirg by corresponding baryon spectral 
functions. B = A^(940), A(1232) are hr upper panel whereas 
B = A*(1620), A*(1650),iV*(1720) are in lower panel for 
fixed values of three momentum of A (fc = 0), temperature 
(T = 0) and baryon chemical potential (^jv = 0). 



FIG. 3: Same as Fig. ([21 for the rest of the baryons 
B = A*(1440), A*(1520), A*(1600) (upper panel) and B = 
A*(1535), A*(1700) (lower panel). 


ending of the dashed line, the Landau regions for differ¬ 
ent loops are clearly visible. As an example for ttTV loop 
the Landau region is M = 0 to toat — i.e., 0 to 0.8 

GeV. Due to the folding of the baryon spectral functions, 
these sharp endings are smeared towards higher value of 
M. Since S'^(M) also depends on T, and k therefore 
total contribution of S'^(M) from all the loops has been 
shown in Fig. (|1]) for different sets of T, /tat and k. 

The nucleon thermal width Fat is basically the con¬ 
tribution of ImE^ at M = ruAT, which is marked by 
dotted line. Being an on-shell quantity. Fat is associ¬ 
ated with the thermodynamical probability of different 



FIG. 4: Imaginary part of total self-energy for different sets of 
nucleon momentum (fc), temperature (T) and baryon chemi¬ 
cal potential (pjv). 



FIG. 5: The variation of nucleon thermal width Fat (upper 
panel) and its corresponding mean free path Ajv (lower panel) 
with T are shown. 


on-shell scattering processes instead of off-shell scatter¬ 
ing processes as described W Weldon for the imaginary 
part of self-energy in Ref. [^. Following Weldon’s pre¬ 
scription, forward and inverse scattering of nucleon can 
be respectively described as follows. During propagation 
of TV, it can disappear by absorbing a thermalized tt from 
the medium to create a thermalized B. Again N can ap¬ 
pear by absorbing a thermalized B from the medium as 
well as by emitting a thermalized tt. The n;(l — n)[) and 
n+(l -bn/) are the corresponding statistical probabilities 
of the forward and inverse scattering respectively [s^, 
because just by adding them, we will get the thermal 
distribution part of Eq. (l28ll . i.e., (n/ +n^). 

From Eq. ((Ml) or (IMl) , we see that Fat depends on 
temperature T, baryon chemical potential /^Af and three 
momentum k of nucleon. The upper panels of Fig. m 































7 



FIG. 6: The variation of nucleon thermal width Fjv (upper 
panel) and its corresponding mean free path Ajv (lower panel) 
as a functions of baryon chemical potential are shown. 



FIG. 7: (color on-line) The T dependence of tjn (upper panel), 
sjv (middle panel) and tin/sn (lower panel) of the nucleonic 
component. The straight red line denotes the KSS bound. 


and dni) are, respectively, displaying the variation of F^v 
with T for different sets of (fc, /i^v) and of F^v with 
for different set of {k,T). The mean free path can be 
defined as XN{k,T, ^ -- and its corre- 

rN{k,T,fiN) 

spending variation with T and fiN are respectively shown 
in the lower panels of Fig. m and dH). The range of 
T and fiN, in which Aat is smaller than the dimension 
of the medium (^ 10 — 40 fm, a typical dimension of 
strongly interacting matter, produced in the laboratories 
of Flic), plays the main role of dissipation via scattering 
in the medium because the larger A is associated with 
the scenario after freeze out of the medium. From the 
dashed line of Fig. ([5|) we see that T > 0.120 GeV (but 
up to Tc ~ 0.175 GeV) is that relevant region for baryon 
free nuclear matter (/i^r =0). Whereas for finite baryon 
chemical potential (e.g. solid line of Fig. at /iAT = 0.7 
GeV), this relevant T region will be shifted slightly to¬ 
ward lower temperature (in addition, Tc is also expected 
to decrease with increase of /tat). Since high momentum 
(fc) of constituent particles always helps them to freeze 
out from the medium, the relevant T region for nucleon 
with high fc is reduced by shifting towards the high T 
region. This can be understood by comparing the solid 
and dotted lines in the lower panel of Fig. ®. 

Using the numerical function T^{k,T, in Eq. (I23L 
we get 77 at as a function of T and which are shown in 
the upper panels of Fig. © and ®. Here we see is 
monotonically increasing with T and both. Using the 
simple equilibrium expression of entropy density (sn) for 
nucleons, 

SN =4:13 J - fXN^n+{uj^), (35) 

the tjn/sn has been generated as a function of T and 
fXN. From the lower panels of Fig. © and m , we see 



FIG. 8: (color on line) The variation of t^at (upper panel), 
Sjv (middle panel) and »7jv/sjv (lower panel) of the nucleonic 
component with pjv. 



FIG. 9: The contributions of different ttB loops in r;jv(r) (left 
panel) and T7jv(/rjv) (right panel). 






























































0.0015 



T (GeV) 

FIG. 10; The T dependence of shear viscosity for pionic (dot¬ 
ted line), nucleonic (dashed line) components and their total 
at piv = 0 (solid line) and /rjv = 0.3 GeV (dash-dotted line). 
The upper and lower panel contain the results without and 
with mixing effect, obtained from Eq. and 1391) respec¬ 
tively. 



FIG. 11: Corresponding results of Fig. (nni) against the /rjv 
axis with two different temperatures. 


that tjn/sn can be reduced by increasing T as well as by 
decreasing /tm. 

In the left and right panels of Fig. the contribu¬ 
tions of different loops (dominating loops only) are indi¬ 
vidually shown in yyjv vs T and ijn vs /i^v graphs respec¬ 
tively. The ttA loop plays a leading role to generate the 
typical values (0.0001 — 0.01 GeV^) of rjN for strongly 
interacting matter because the major part of the nucleon 
thermal width is coming from this loop only. 

Up to now, we have calculated the contribution of shear 
viscosity from nucleon thermal width, although a ma¬ 
jor contribution comes from the thermal width of pion. 
Hence, one should add the pionic contribution with nu¬ 
cleon contribution for getting total shear viscosity of nu¬ 
clear matter at finite temperature and density. In our 


recent work [^ . the shear viscosity, coming from pionic 
thermal width has already been addressed. The one-loop 
Kubo expression of shear viscosity and ideal expression 
of entropy density for pionic components are respectively 
given below. 



where nfc(w^) = l/{e^“^ — 1} is the Bose-Einstein distri¬ 
bution function of pion with wj = (P -|- and F^ 

is the thermal width of tt mesons in the medium due to 
TTCT and Tip fluctuations. 

Now, adding that pion contribution with the nucleon, 
one can simply get the total shear viscosity of nuclear 
matter as 


Vtot — ^7r T VN ; 


(38) 


where and do not face any mixing effect of pion 
density, = 3 f ■^|^nfc(w^) and nucleon density, pjv = 

However, viscosity of single component 
gas should be different from the viscosity of that compo¬ 
nent in a mixed gas [s^ l . To incorporate this mixing 
effect for rough estimation, we follow the approximated 
relation (^. [54| 

= (39) 


where 


_(Ft_ 

i + (^) (?e)\P¥™ 


and 


Vn 


_ m _ 

1 _L f P2L.\ ( \ , / l-t-mw/mT 

^ + Vpn J ) y 2 


(40) 


(41) 


For simplicity, the cross sections of all kinds of scattering 
are taken as constant with same order of magnitude {i.e. 

~ cf-kN ~ o-jvAf). In the upper panels of Fig. CQD 
and (ED, the T and pN dependence of (dotted line), 
Pn (dashed line) and their total ptot (solid line and dash- 
dotted line for two different values of pN and T) are sep¬ 
arately shown. Whereas lower panel of the figures show 
their corresponding mixing effect following from Eq. (HOI) . 
(ITTll and (1551) . Erom the Fig. (ITT]) , one should notice that 
the independent nature of p-^ipN) has been changed to 
a decreasing function due to mixin g eff ect. Similar qual¬ 
itative trend has been seen in Ref. 241. 

The entropy density of nucleon component from 
Eq. (13511 . pion component from Eq. (ED and their to¬ 
tal Stot = SAf + s,r are individually shown in the up¬ 
per panels of Eig. (ED and ED as functions of T and 
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FIG. 12: (color on-line) The entropy density (upper panel), 
viscosity to entropy density ratio withont (middle panel) and 
with (lower panel) mixing effect as functions of T. 



FIG. 13: (color on-line) The corresponding results of Fig. (1121) 
are shown with respect to /tat. 


fiN respectively. The corresponding 77 /s without (middle 
panel) and with (lower panel) mixing effect are shown 
in Fig. (HU and m as a function of T and hn respec¬ 
tively. The decreasing nature of total rj/siT) qualita¬ 
tively remains the same after incorporating the mixing 
effect whereas an increasing function of the total 77 / 5(71 at) 
transforms to a decreasing function due to this mixing ef¬ 
fect. Comparing our results with the results of Itakura 
et al. [l^, where rj/ s{^n) also reduces with the mix¬ 
ing effect appears to be very important. However, the 
total ri/s{fiN) in mixing scenario becomes an increasing 


function beyond /ijv ~ 0.5 GeV because the increasing 
rate of dominates over the decreasing rate of 

in that region. Using the effective hadronic 
Lagrangian, the conclusion of our results should be con¬ 
centrated within regions of 0.100 GeV <T < 0.160 GeV 
and 0 < < 0.500 GeV. 

V. SUMMARY AND CONCLUSION 

Owing to the Kubo relation, the shear viscosity can 
be expressed in terms of two point function of the vis¬ 
cous stress tensors at finite temperature. By using the 
real-time thermal field theoretical method, this two point 
function has been represented as NN loop diagram when 
the nucleons are considered as constituent particles of the 
medium. A finite nucleon thermal width Fat has been 
traditionally included in the nucleon propagators of the 
NN loop for getting a non-divergent shear viscosity 77 at. 
This nucleon thermal width is obtained from the one-loop 
self-energy of nucleon at finite temperature and density. 
Different possible pion baryon loops are accounted to cal¬ 
culate the total Fat, which depends on the three momen¬ 
tum of nucleons (k) as well as the medium parameters 
T and Using the numerical function rAr(fc, T, TiAr), 
77 AT and tjn/sn are numerically generated as functions of 
T and un- Adding the pionic contribution taken from 
Ref. |53l | with the numerical values of the nucleonic com¬ 
ponent, we have obtained the total shear viscosity, where 
a gross mixing effect of two component system has been 
implemented. Along the temperature axis, the shear vis¬ 
cosity of both pion and nucleon components appear as 
increasing function, whereas along the /^Af axis shear vis¬ 
cosity of pion component changes from its constant be¬ 
havior to a decreasing function due to presence of mixing 
effect. The total shear viscosity to entropy density ratio 
(^tor/stot) for the pion-nucleon mixed gas reduces with 
increasing T as well as and quantitatively becomes 
very close to the KSS bound. This behavior indicates 
that 77j((},’'/stot tends to reach its minimum value near the 
transition temperature at vanishing as well as finite value 
of /7Af- According to these results, the finite baryon chem¬ 
ical potential helps the nuclear matter to come closer to 
its (nearly) perfect fluid nature. 
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